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When radio rays pass through the atmosphere, they are defocused due to its presence. 
This effect is measured by the divergence coefficient and general formulas are derived for Di, 
the divergence coefficient of the direct ray, and Z>2, the divergence coefficient of the reflected 
ray — assuming a smooth spherical earth. As examples, Di and D2 are shown for some 
typical cases with an ''exponential" atmosphere (troposphere). 

1. Introduction 

Tlio field at a radio receiver within line-of-sight of the transmitter may be found by adding 
vectorially tlie field due to the direct ray and tliat due to the reflected ray. These fields depend 
directly upon the defocusing of the respective rays and this effect is measured by Z^i, the 
divergence coefficient of the direct ray, and A, the divergence coefficient of the reflected ray. 

The direct ray is tlie ra}^ which goes directly from the transmitter to the receiver and it 
is defocused due to the change of refractive index, n, with height. This defocusing effect is 
given by : 

A=|-; (1.1) 

where Ei is the strength of tlie electric field due to the ''direct'^ ray at the receiving point, P (in 
the atmosphere), and E'o is the field strength which would be observed at the same distance if 
the system were located in a medium of constant refractivity. 

The reflected ray strikes the earth (assumed spherical) and is reflected to the receiver. 
It receives some additional defocusing because of the earth's shape, and the total effect is 
given by : 

A=§ (1.2) 

where E2 is the strength of the electric field due to this ray at the point P and E is the field 
strength that would be observed at that distance if the system were in a homogeneous medium 
and the earth were flat. 

2. Divergence Coefficient of the Direct Ray 

The divergence coefficient of the direct ray, Di, may be derived by first considering a radio 
ray passing through the atmosphere unimpeded. The energy of this ray is 

Energy = Cn2E\dq, = Cn.Eldqo (2 . 1) 

where C is a constant of proportionality, Ei is the field strength at the receiving point P — • 
where the index of refraction is /I2 — and dqi is the cross section of the ray at this point. Eq and 
dg[o are the corresponding quantities at the same distance if the system is in a medium whose 
index of refraction is Ui, that at the transmitter. 
Then from (1.1) 
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Figure 1. Geometry of the direct ray. 



Referring to figure 1 



dqo= 



and 



tRKcIM^ 



r/?o cos I3id^id<t) 



dqi= 



TTpl sin ^2 sin doddodcf) 



(2.3) 
(2.4) 



where is the angle the ray makes with the plane of the paper. 
Combining (2.2), (2.3), and (2.4): 



P2 \ I 



rii cos jSi 



712 sin 132 sin do{ddo/dl3i) 



(2.5) 



^0 is the central angle and, for a radio ray passing through a spherically stratified atmos- 
phere, it is given by [Kerr, 1954; Counter, 1956]: 



where the constant K is 



'•=1 f' 



dp 



K=ntPi cos |8j, 



(2.6) 



(2.7) 



which is Snell's Law. (See also appendix I.) In the above equations, p is the distance from 
the center of the earth to a point on the ray and n is the index of refraction at that point; Pi and 
Ui are the corresponding quantities at the end points and ^t are the angles (— 7r/2</3i<7r/2) 
between the ray at these points and the horizontal. In this paper, t = l is associated with the 
transmitter and t = 2 with the point P, although they may be interchanged in the final result, 
(2.5). 

The derivative of (2.6) with respect to j3i gives 



\deo/d0,\ = 



tan 



K' 






n^pdp 



[{np/Kr-l] 



3/2 



(2.8) 



(Again see appendix I.) 
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The distance 7?o is not completely defined; e.g., it may be the distance along the ray or the 
straight line distance to the point P — the only requirement being that 



/io= 



^[p?+pl— 2piP2 COS 6^ 



(2.9) 



when n is constant, since the ray becomes straight. In the calculations described below, i?o 
is the straight line distance and therefore is defined by (2.9) for variable n also. 

For constant n the integrals in (2.6) and (2.8) may be evaluated in closed form and, putting 
the results in (2.5), A^^l, the value which is usually used. 

In the troposphere n is of the form [Bean and Thayer, 1959] 



n=--l + {N,Xl()-')e- 



(2.10) 



where Ns is the refractivity on the earth's surface, c is the decay constant, and h=p—po is the 
height above the earth; po is the radial distance to the earth's surface and the following empirical 
relation reflects the average correspondence between Ns and po : 



Po (km)-6370+10ln [1 + 0.6 exp (-3.35X10-^WJ)] 



(2.11) 



where 6,370 km is taken as the earth's radius. The refractivity decay constant c in (2.10) 
is given by [Bean and Thayer, 1959; Rice, Longley, and Norton, 1962]: 



c (per km) = H 



^"[iv:=7:3 



iV« 



J 



;iV>250 



.32 exp (0.005577iV,) 
A^,X10-*[7.939-0.01166iV,];iV,<250 



(2.12) 



Table 1 lists values of po for Ns=200, 300, and 400. Using these values, the integrals in (2.6) 
and (2.8) were evaluated as sliown in appendix II, Eq was computed from (2.9), and Di was 
plotted logarithmically against the takeoft' angle ^i in figures 3 and 4. In figure 4, the cutoff is 
due to interference by the earth. 

In these figures, hi is the height of the transmitter above the earth's surface (0 and 5 km) 
and fh is the height of the receiving point, P, above the earth's surface: 



h—Pi — po', 

h2=p2—PO' 



(2.L3a) 
(2.13b) 



Table 1 



Ns 


Po (km) 


c (per km) 


200 
300 
400 


6, 373. 008823 
6. 370. 390117 
6, 370. 001131 


0. 11214 

. 1392842847 
. 1867197187 



3. Divergence Coefficient of the Reflected Ray 



The divergence coefficient of the reflected ray, D2, may be derived by again considering a 
radio ray passing through the atmospliere but now reflecting from the earth before reaching 
the receiving point, P. Kei'erring to figure 2, Do is. developed in the same manner as Di and the 
result is 



P2 V 



D,= 



Til COS /3i 



ri.2 sin J82 sin d{dd/dfii) 
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(3.1) 




Figure 2. Geometry of the reflected ray. 



The central angle B is now the sum of the angles di and 62] 

where ^1 and 62 are [Kerr, 1954; Counter, 1956]: 

^^IPi dp 

and 

2 IJpij „f/'«„/S^2_lll/2 



^2 



dp 



;Po p[(np/i^^-l] 
if the atmosphere is spherically stratified, if is the same as before, i.e., 

K=niPi cos ft 



(3.2) 

(3.3a) 
(3.3b) 

(3.4) 



where pt is again the distance from the center of the earth to a point on the ray, rii is the index 
of refraction at that point, and Pi is the angle the ray makes with the horizontal at that point; 
i=0 is associated with the earth's surface, i=l again with the transmitter, and i=2 with the 
point P, although these may be interchanged in (3.1). 
The derivative of 6 with respect to /5i then gives 



\d6/dp^\ = 



tanj3 



" \Jpo [(:np/K')-ir'^J.o [{np/K)'-ir'f[ ^^'^^ 



R is defined only by the fact that it must be the slant range, 

R=\pl+pl-2pop, cos d,Y^'+[pl+pl-^PoP2 cos 62]'^', 
when n is constant. 



(3.6) 
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For the above case {n constant), the curvature of the earth is the only contributing factor 
in determining D2, and 



' P2\sin(9(i 



sin /3o cos^ jSi 



sin ^2 sin ^i+sin ft sin 6^ 



(3.7) 



which is, as expected, the same as that derived by Van der Pol andlBrenmier [1939] and Riblet 
and Barker [1948]. 

For small heights, pi and P2^Po- (3.7) may be approximated by 



26,62 1"^ 



^^\_ "^eitaniSoJ 



(3.8) 



This is the formula which is usually used with the earth replaced by one with a radius 
k{pQ) times the actual radius [Norton, 1941], where 



^(Po) ' 



n 



~n+p{dn/dp) 



(3.9) 



It is evaluated for some typical cases with this assumption and the results are shown in figure 5, 
together with values of (3.7) without this assumption; in this figure, results are also shown for 
the same situations but assuming an exponential atmosphere (troposphere) with A^s=200, 300, 
or 400. (See (2.10) and table 1.) The transmitter height hi is 5 km. 

It should also be noted that the curves in this figure are plotted against the reflection angle 
l3o instead of the takeoff angle ft. ft may be found, however, by using Snell's Law, (3.4), and 
ft, for a particular path, may be found by applying an iterative method to (3.2). 

4. Conclusion 

The defocusing effect of the atmosphere may be of importance for small takeoff angles, 
ft. (See figs. 3 and 4.) This defocusing effect, however, is usually small compared to that 
arising from the earth's curvature. This is indicated by comparing the values of D2 (atmos- 
phere plus the earth's curvature, fig. 5) with the corresponding values of Di (atmosphere only, 
fig. 3). 

Especially for large antenna heights, a good approximation of A, the divergence coefficient 
of the reflected ray, is then given by (3.7)^no atmospheric effect — with ft being the reflection 
angle of the true ray. This is also shown in figure 5. 
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Figure 3. The divergence coefficient of the direct ray with one terminal on the ground. 
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Figure 4. The divergence coefficient of the direct ray with one terminal elevated S k7n. 
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Figure 5. The divergence coefficient of the reflected ray with one terminal elevated 5 km. 

5. Appendix I 

When i3i<0, there may be two points on a ray at the same height and, for iS2>0, (2.6) 
should be replaced by the sum 






dp 



plinplKY-lf'^ ' J, p[(np/i^)2-l]>^ 
where r is such that rn{r)=K. Then (2.8) should be replaced by 



\ddom\ = 



J . . cos ^2 sm ft w N , . ^ f^i {dkdp)dp . ^ ^ C 

Hpi) ^—- — o ^(P2) + tan ft r/ /r^xV .11/ + tan ft 

cos ft sm ft ^"^^^ ' "^'J, ?2p/i^)2— 1 3/2 ' ^ij 



cos ft sin ft "^^^^^ ' ^ ^'J, [(r2,p/i^)2-l] 
where k(pi) is again defined as [CCIR Study Group IV-C, 1959]: 

HPi)^: 



P2 {dk/dp)dp 
[{np/Ky-l] 



n+p{dnjdp) 
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(I.l) 



(1.2) 



(1.3) 



6. Appendix II 

In finding do and \ddo/d^i\ for an exponential atmosphere, the integrals were evaluated by 
first letting r cos z = pi cos ^i and then using Gaussian quadrature. 

The only times when this method could not be used were in finding \ddoldP\\ wlien i3i = 
or ^2 = ^' For these cases integration by parts was employed, with np being tlie variable, and 
it was found that 

\ddjd^,\=k{p,) (ll.l) 

when j(3]=0, and that 

|sin ^2{de,ld^,)\ = \k{p2) tan /3i| (11.2) 

when j(32 = 0. For the definition of k{p^) see appendix I. 

7. References 

Bean, B. R., and G. D. Thayer, On models of the atmospheric radio refractive index, Proc. IKE 47, No. 5, 

740-755 (May 1959). 
CCIR Study Group IV-C, CCIR Athises of propagation curves, resohitions 21 and 22 (Geneva, 1959). 
Counter, V. A., Propagation of radio waves through the troposphere and ionosphere, Lockheed Aircraft Co., 

Missile Systems Div., Palo Alto, Calif. (Dec. 1956). 
Kerr, D. E., Propagation of short radio waves, p. 49 (McGraw-Hill Book Co., Inc., New "^'ork, N.Y., 1954). 
Norton, K. A., The calculation of ground-wave field intensity o\'er a finitely conducting s[)hcrical earth, Proc. 

IRE39, 623 (1941). 
Riblet, H. J. and C. B. Barker, A general divergence formula, J. Appl. Phys. 19, 03-70 (Jan. 1948). 
Rice, P. L., A. G. Longley, and K. A. Norton, Transmission loss prediction for tr()})ospheric communication 

circuits, to be published as an NBS Technical Note in 1962. 
Van der Pol, B., and H. Bremmer, Further note on the propagation of radio waves over a finitely conducting 

spherical earth, Phil. Mag. 27, sec. 7, No. 182 (Mar. 1939). 

(Paper 66D4-209) 



485 



